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ABSTRACT 

We present the results of a kinematic study of planetary nebulae in the extreme outskirts of two 
spiral galaxies, M83 (NGC 5236) and M94 (NGC 4736). We find that in the inner regions of the 
galaxies, the vertical velocity dispersion (a z ) falls off exponentially with the light, as expected for a 
constant mass-to-light ratio, constant thickness disk. However, starting at four optical scale lengths, 
cr 2 asymptotes out at roughly 20 km s _1 . Our analysis finds evidence for significant flaring in the 
outer regions as well, especially in M94. These observations are in excellent agreement with predictions 
derived from models of disk heating by halo substructure, and demonstrate how kinematic surveys in 
the outer disks of spirals can be used to test hierarchical models of galaxy formation. 
Subject headings: galaxies: individual (NGC 5236, NGC 4736) — galaxies: kinematics and dynamics 
— galaxies: spirals — planetary nebula: general 



1. INTRODUCTION 

The cold dark matter (CDM) paradigm has proven 
to be very successful in explaining the la rge-scale struc- 
ture of the universe and galactic clusters (Tcgmark ct al. 
I2004T ) . yet problems still exist at the galactic level. One 
such issue is the "missing satellite" problem: according 
to numerical simulations, the structure of a galactic halo 
should look like a small version of a galaxy cluster and 
contain many sma ll subhalos (e.g.. [Moore et ail 119991 : 
iKlypin et al.| [l999). The recent discoveries of satellite 
galax i es tidally stripped by the Milky Wav (pbata et all 
119941: iNewberg et ail l2002t Uuric et all 120081. ari d reL 
erences therein) and Andromeda (|Ferguson et ail 120051 : 
llbata et aLll2007f ) help to solve this problem. These small 
companions have a profound effect on the morphology 
and kinematics of thin galactic disks. Specifically, nu- 
merical models of satellite halo bombardments predict 

1. the formation of faint stellar streams above the disk 
plane (which would be very difficult to observe in 
any but the nearest galaxies), 

2. the generation of long-lived, low-surface brightness, 
dynamically cold, ring-like features in outer disks, 

3. the growth of a strong bar, 

4. the production of a pronounced flare, and 

5. the development of a thick disk 



(jKazantzidis et ail 120081 . and references therein). 
Though many simulations (iQuinn e~ al. 1993; Fon t et ail 
120011 llavashi fc Chibal [20061 : iKazantzidis et all [20081 . 

and references therein) have been performed to explore 
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disk heating by halo substructure, observationally the 
results from this process have proven to be rather elusive. 
Only recently have teams begun to d etect signs of this 
proce ss in the Andromeda Gal axy ([Richardson et ail 
l2008f) and other nearby galaxies (|de Jong et al.ll2008f ). 

Here we describe kinematic evidence for halo substruc- 
ture derived from the velocity dispersion of planetary 
nebulae (PNe) in the extreme outer disks of two nearby, 
face-on spirals. We use stability arguments to show that 
the disks of M83 (NGC 5236) and M94 (NGC 4736) must 
flare at radii greater than ~4 scale lengths and demon- 
strate that at these large radii, the z velocity dispersions 
of the old disk stars agree with the results of numer- 
ical simulations of a di sk heated by halo substructure 
(|Havashi fc Chibal [200l . 

2. THE SURVEY 

Flat rotation curves indicate the presence of dark mat- 
ter in the outer regions of spiral galaxi es and allow us 
to de termine total galactic mass (e.g., ISofue fc Rubinl 
2001). However, rotation curves alone cannot decou- 
ple the m ass contribution of the disk from that of the 
dark halo (jBarnes et al.ll2004l ). To break this degeneracy, 
we have been using the motions of planetary nebulae in 
low-inclination spirals to measure disk mass directly via 
the z motions of stars. PNe are ideal particles for this 
purpose: they are bright, abundant to >5 scale lengths, 
representative of the old disk, relatively easy to distin- 
guish from H II regions (via their distinctive [O III]-Ha 
ratio; Ciardullo et al. 2002), and amenable to precise 
(~2 km s _1 ) radial velocity measurements with fiber-fed 
spec trographs. 

In Herrmann ct al. (2008) (Paper I), we presented the 
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results or a narrow- band imaging survey of six low- 
inclination nearby spirals (IC 342, M74 (NGC 628), 
M83 (NGC 5236), M94 (NGC 4736), NGC 5068, and 
NGC 6946) in which we identified 165, 153, 241, 150, 
19, and 71 PNe, respectively. In two upcoming papers 
(K.A. Herrmann & R. Ciardullo 2009a, 2009b, in prepa- 
ration, Papers II and III), we will detail our spectroscopic 
follow-up observations, present high-precision radial ve- 
locities for 550 of our PN candidates, and use the data 
to estimate dynamical disk masses. 



2 



Herrmann, Ciardullo, & Sigurdsson 



Here, we focus on the data set for the two galaxies for 
which we have the largest radial coverage, the SBc spiral 
M83, which is at a distance of 4.8 Mpc, and the earlier 
Sab system M94 at D = 4.4 Mpc (Paper I). Our PN ve- 
locity sample for the former galaxy consists of 162 objects 
at radii between 2.5 and 24.7 kpc, i.e., between ~1 and 10 
optical disk scale lengths. These data, which were taken 
with the Hydra bench spectrograph on the Cerro Tololo 
Inter-American Observatory (CTIO) 4m telescope, have 
typical velocity uncertainties of a v ~ 6.5 km s _1 , and in 
all cases, a v < 15 km s _1 . For M94, our data set contains 
127 planetaries observed with the Hydra spectrograph 
on the WIYN (Wisconsin, Indiana, Yale & NOAO) tele- 
scope. These PNe have galactocentric radii between 0.5 
and 8.8 kpc (0.4 < R < 7 disk scale lengths), and typi- 
cal velocity uncertainties of a v ~ 4 km s , again with 
a v < 15 km s^ 1 for all objects (see Paper II). 

Since neither galaxy is exactly face-on, we began our 
analysis by removing the effects of galactic rotation from 
the PN sample. This was done using veloci ty maps from 
The H I Nearby Galaxy Survey (THINGS; IWalter et all 
|2008|) with a correction for asymmetric drift. We then 
binned the planetaries by radius (with 15 — 16 PNe per 
bin in M83, and 17 - 18 PNe per bin in M94), and iden- 
tified those objects more than ~2.5 a away from their 
bin mean as possible halo contaminants. (In practice, 
this made very little difference to the analysis, since the 
procedure eliminated only six objects in M83 and three 
in M94.) Finally, to extract the component of the ve- 
locity dispersion perpendicular to the galactic disk, a z , 
from the other two constituents of the velocity ellipsoid, 
we began by using the epicyclic approximation for near- 
circular orbits to couple a$ to <tr. We then constrained 
the shape of the disk velocity ellipsoid using the lim- 
its imposed by the physics of disk scatter ing {a z < or; 
IVillumser] 119851 : I Jenkins fc Binnevl 1 19901) and bending 
instabilities {a z > 0-25(7^; Hbomri U966j lArakil 119851 : 
iMerritt fc Sellwoodlll994D , and computed the likelihood 
that each combination of a z and an could produce the 
line-of-sight dispersions observed within our bins (see Pa- 
per III). 

Figure [T] shows the results of this analysis, display- 
ing contours which enclose 38% (0.5<r), 68% {la), 86% 
(1.5cr), and 95% {2a) of the probability. Note that these 
contours have a slight tilt to them. This is a consequence 
of the galactic inclination: neither galaxy is precisely 
face-on {i ~ 24° for M83, i ~ 35° for M94), and the larger 
the inclination, the more difficult it is for our maximum 
likelihood procedure to extract a z from the line-of-sight 
velocity dispersion. Nevertheless, the figure does demon- 
strate that our measurements of a z are reasonably well 
defined, especially for M83. 

3. RESULTS 

For spiral galaxies, the velocity dispersion perpendicu- 
lar to the disk is related to disk surface mass, E(-R), and 
the disk scale height, h z , by 

a 2 z {R) = KGT,{R)h z , (1) 

where K = 2ir for the isothermal case, K = tt 
for disks whose vertical mass density drops exponen- 
tially, and K = tt 2 / 2 for the intermediate sech(2:) case 
van der Kruitlfl98l . If (1) the disk mass-to-light ra- 
tio is constant, (2) the disk light of the galaxy decays 
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Fig. 1. — Maximum likelihood probability contours for the outer 
disks of M83 and M94. From dark to light, the contour regions 
enclose 38% (0.5cr), 68% (la), 86% (1.5fr), and 95% (2cr) of the 
probability. The dashed lines show the limits of our analysis; the 
black crosses indicate the most likely solutions. The solid curves 
display the upper limits on <r z derived from the additional con- 
straint that the disk be stable against axisymmetric perturbations. 
The multiple curves represent intermediate (sech (2)) disks with 
scale heights of 200 and 300 pc. Note that for stability, h z must 
be larger than these nominal values, especially in M94. 

exponentially with a single scale length {ha), and (3) 
the disk scale height is also constant, then a z should de- 
crease exponentially, with a scale length twice that of 
the light. As the marginalized probabilities of Figure [2] 
illustrate, this is not the case in either galaxy. While 
a z does track the light in the galaxies' inner regions, 
the curve flattens out at distances more than four disk 
scale lengths from the nucleus. The minimum disper- 
sion of a z ~ 20 km s _1 is much greater than the typical 
measurement error of <5 km s , and greater than the 
<15 km s" 1 values expected from t he velocity dispersion 
of ionized gas an d H II regions (i.e.. lZaritskv et al]|1990l : 
iFathi et al.ll2007l ). Moreover, it is difficult to conceive of 
any way that internal extinction could create such an ef- 
fect. While it is true that nonisothermal disks may have 
velocity di spersions that incre ase with height above the 
plane (see Ivan der Krurnll988h . any layer of dust that is 
thick enough to effect our derived values of a z would also 
extinct the bulk of the PN population below the detec- 
tion threshold and alter the distribution of PN [O III] /Ha 
line ratios (see Paper II). Finally, large-scale warping of 
the disks is not a solution: even with modeling isophotal 
twists by adjusting the disk position angle and inclina- 
tion by ~70° and ~15°, respectively, we cannot produce 
the dispersion curve seen in Figure O (see Paper III) . 

The flat velocity dispersion profile of M94 is partially 
explained by its unique 3.6 /im-band radial profile, which 
has one scale length {Hr = 1.22 kpc) in its inner regions, 
then breaks to follow a shall ower profile {Hr = 7 .16 kpc) 
at radii greater than 5 kpc (|de Blok et al.|[2008h . (A re- 
cent deep i?-band profile shows the same broken shape 
(Erwin 2009, private communication).) But M83's disk 
has no such break, as its simple exponential profile 
{fiR = 2.45 kpc) extends over ~20 k pc to the edge of 
our survey area (|de Jong et al.l f2008h . Moreover, even 
with its broken disk profile, M94's a z values are much 
higher than possible for a constant M/L, constant h z 
disk. 

These high a z values also bring up another issue. To be 
stable against axisymmetric perturbations, a thin stellar 
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disk must obey the iToomrd (|1964f ) criterion 



where k is the epicyclic frequency of the orbits. Combin- 
ing this criterion with Equation {1} yields a constraint 
on the allowable values of a z : 

( Kh z K<r B .\ 1/2 , . 

Because our galaxies are low-inclination systems, h z is 
not directly obtainable from our observations. However, 
imaging surveys of large, edge-on, noninteracting spi- 
rals demonstrate that the scale length to scale height 
ratio of galactic disks ranges from ~10 in late-type sys- 
tems to ~5 i n ear lier-type objects (e.g., Ide Griisl Il998t 
iKregel et al.l 120021 ). When applied to the two galax- 
ies studied here, this relation implies that h z should be 
between 200 and 300 pc for both t he Sab spiral M94 
(h R = 1.22 kpc; Ide Blok et ail 12008ft and the SBc sys- 
tem M83 (h R = 2.45 kpc; Ide Jong et al.ll2008ft . Yet, as 
the curves of Figure Q] indicate, outer disks as thin as 
this are unlikely, as our solutions for a z in M83 and es- 
pecially M94 lie predominantly above the stability limit. 
This suggests that the stellar disks of these systems flare 
dramatically in their outer regions. 

Could this flaring and higher than expected values of 
a z be due to heating of the disk by halo substructure? 
We can address this question by comparing our results 
to iV-body simulations of dark subhalo interactio n s with 
large galaxies. In particular, [Hayashi fc Chibal (2006) 
have modeled the heating of a Milky Way-like galaxy 
b y a population o f subhalos, distributed according to 
a iHernquistl ([1990) law. Figure [2] compares our a z ve- 
locity dispersion profile to their Model F, which repro- 
duces interactions between an initially thin, stable, con- 
stant mass-to-light ratio disk, and a few (~300) massive 
(<10 9 M Q ) subhalos, totaling M sh < 0.l5M disk and dis- 
tributed with a half-mass radius of ^210 kpc. Although 
their model only probes out to ^5 disk scale lengths, 
the agreement between observations and theory is ex- 
cellent. Their model not only reproduces the run of a z 
versus radius seen in our galaxies, but also predicts an 
amount of disk flaring that is c onsist ent with our anal- 
ysis. Moreover, IStrigari et all (|2008ft have recently re- 
ported that, of the 23 known Milky Way satellite galax- 
ies, the 18 with dynamical mass measurements all have 
masses of M ~ 10 7 M Q within a radius of 0.3 kpc. This 
suggests that these satellites once had total masses of 
^10 9 M P) before enco untering the Milky Way's potential 
(Stri gari et alJl2008f ). This value is i n excellent agree- 
ment with the subhalo masses used in lHavashi fc Chiba 
(|2006ft . and supports the idea that interactions with sub- 
halos could be responsible for the flat dispersion profiles 
seen in the outskirts of galaxies. 

4. CONCLUSIONS 

The existence of persistent very thin disks in spiral 
galaxies has long been recognized to be a si gnificant con- 
strain t for models of structure growth (cf. iMoore et al.l 
[19991) . Major mergers destroy thin disks, while large 
asymmetries or substructure in extended dark matter 
halos will warp or heat these structures (|Quinn et al.l 



119931: iFont et all 1200 ll ). Conversely, evidence for ex- 
tended thick disks (or thin young disks embedded within 
older, thicker disks) can preserve direct evidence of past 
mergers and/or episodes of externally triggered star for- 
mation. 

The two galaxies considered here have very different 
radial profiles: while the slope of M94's "antitruncated" 
stellar disk changes dramatically at a radius of ~5 kpc 
(|de Blok et al.ll2008| ). M83's exponential profile appears 
undis turbed out to ~10 disk scale lengths (|de Jong et al.l 
2008). Thus, while M94 may h ave been involved i n 
a minor merger some years ago ([Younger et al.ll2007t ). 
the only evidence for interaction in M83 is the warp- 
ing of its extreme outer H I disk, which begins ~8 scale 
lengt hs from the nucleus (|Huchtmeier fc Bolmen stengell 
Il981h . Yet in both these systems, the kinematic evi- 
dence suggests that the stellar disk flares dramatically 
past ~4 disk scale lengths, and that at these large dis- 
tances, the z velocity dispersion is independent of radius. 
This result is consistent with cosmological models of hier- 
archical structure formation, where a thin disk is heated 
by a relatively small number of massive subhalos, which 
are either embedded within the parent dark halo with a 
relati vely small initial spatial spre ad (|Havashi fc Chibal 
12006ft . or initially on radial orbits ([Hopkins et al.l l2008). 
Since simulations indicate the disk heating does not scale 
linearly with the mass and number of subhalos, this is 
potentially a strong constraint on substructure. While 
there are several different processes that may combine to 
heat a disk, cooling such a disk back to a low velocity 
dispersion state is difficult. (The slow accretion of ex- 
tremely large amounts of cold gas could, in theory, lead 
to adiabatic compression of the disk. However, in such a 
process, the gas mass would need to be many times that 
of the stellar disk, and even then, there would be sta- 
bility issues.) Our current measurements of a z at large 
radii already rule out some models for halo substructure, 
and are directly consistent with models where there are 
fewer, more massive subhalos on orbits with relatively 
low angular momentum. 

More data are required to explore this constraint fur- 
ther. In particular, measurements over a wider range 
of galaxy masses, radii, and Hubble types are needed. 
Fortunately, such data are relatively straightforward to 
obtain with moderate-sized telescopes. Planetary nebula 
surveys in the extreme outer disks of the low-inclination, 
late-type spirals M74, M101, and IC 342 are already un- 
der way, and observations in earlier systems, such as M95 
and NGC 1291 are possible. Thus, a robust comparison 
to hierarchical models of structure formation should be 
possible in just a few years. 
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sion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space 
Administration. This work was supported by NSF grant 
AST 06-07416 and a Pennsylvania Space Grant Fellow- 
ship. 
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Fig. 2. — Compariso n of our M83 and M94 results to the a z values (solid black curve) determined from the Model F (hn = 3.5 kpc) 
numerical simulation of Hayashi & Chiba (2006). The data have not been rescaled in any way. The agreement between the data and model 
is much better than with the constant mass-to-light ratio, constant scale- height exponential disk (dashed curve). 
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